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ABSTRACT 

According to the now strongly supported concordance ACDM model, galaxies 
may be grossly described as a luminous component embedded in a dark matter halo. 
The density profile of these mass dominating haloes may be determined by N - body 
simulations which mimic the evolution of the tiny initial density perturbations during 
the process leading to the structures we observe today. Unfortunately, when the effect 
of baryons is taken into account, the situation gets much more complicated due to 
the difficulties in simulating their physics. As a consequence, a definitive prediction of 
how dark matter haloes should presently look like is still missing. We revisit here this 
issue from an observational point of view devoting our attention to dwarf galaxies. 
Being likely dark matter dominated, these systems are ideal candidates to investigate 
the present day halo density profiles and check whether dark matter related quantities 
correlate with the stellar ones or the environment. By fitting a large sample of well 
measured rotation curves, we infer constraints on both halo structural parameters 
(such as the logarithmic slope of the density profile and its concentration) and derived 
quantities (e.g., the mass fraction and the Newtonian acceleration) which could then 
be used to constrain galaxy formation scenarios. Moreover, we investigate whether the 
halo properties correlates with the environment the galaxy lives in thus offering a new 
tool to deepen our understanding of galaxy formation. 
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1 INTRODUCTION 

The discovery of the cosmic speed up through the Hubble 
diagram of Type la Supernovae more than a decade ago has 
pointed at the cosmological constant A as the ideal candi- 
date to drive the accelerated expansion. Almost fifteen years 
later, the ACDM model, made out by the A term dominating 
the energy budget and the cold dark matter (CDM) being re- 
sponsible for the clustering on galaxy and clusters scales, has 
proved to excellently work in fitting a wid e variety of data on 
cosm ological scales, running from SNela ( AmanuUah et al. I 
[2OI0I) to Baryon Acoustic Oscillations (jPercival et al. II2OICI ) 
and cosmic microwave background anisotropy and polariza- 
tion spectra (jKomatsu et al. II2OIH ). 

The unprecedented accuracy in describing the large 
scale universe is frustratingly contrasted by the serious 
shortcomings the ACDM model has to deal with on the 
galaxy scales, the most well known difficulties being rep- 
resented by the substructure problem, i.e. the strong dis- 
crepancy between the predicted and observed number of 



satell ites (|KIvpin et al. 1 1 19991 : [Moore et al. 1 1 19991: iBullock 
201c ) , and the cus p /core controver sy (JFlores fc Primack 
1994 : iMoore 1 1 1994 ide Blok I I2OI0I ') with rotation curves 



data arguing in favour of cored density profiles instead of 
the cusped ones coming out of dissipationless ACDM based 
Nbody simulations. Notwithstanding the significant pro- 
gresses in both modeling and observations during the last 
15 years, both these problems are still on the ground, thus 
motivating further investigations. 

To this end, dwarfs galaxies stand out as ideal candi- 
dates to make our knowledge of haloes a step forward. These 
small mass systems, characterized by having low values of lu- 
minosity, metallicity and size, are important building blocks 
of more massive galaxies. In CDM cosmologies, the smallest 
dwarf galaxies, likely formed before reio nization with masses 
smaller than 10^-10^ Mq (see, e.g., ICnedin fc Ostriken 
Il99'if ). 

Morphologically, dwarf galaxies can be classified as 
dwarf spiral galaxies, blue compact dwarf galaxies (BCDs), 
dwarf irregular galaxies (dirrs), dwarf elliptical galaxies 
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(dEs) , dwarf spheroidal galaxies (dSphs) , ultra - faint dwarf 
spheroidals, and tidal dwarf galajcies. Dwarf spiral galaxies, 
BCDs and dirrs are star - forming objects, gas rich and rota- 
tionally supported, while, on the contrary, late -type dwarf 
spirals are slow rotators or exhibit solid - body rotation. dEs 
usually have little or no detectable gas, and are often not ro- 
tationally supported, similarly to dSphs. Tidal dwarf galax- 
ies, finally, form from the debris torn out of more massive 
galaxies during interactions and mergers. 

Similarly interesting hints on the cusp/core controversy 
may be drawn from Low Surface Brightness (LSB) galax- 
ies. These are late -type, gas -rich, dark matter dominated 
disk galaxies. Their optical appearance is dominated by an 
exponential disk with a young, blue population, with little 
evidence for a dominant old population. Additionally, these 
galaxies do not have large dominant bulgefjj. 

Remarkably, whether the dwarfs haloes are cored 
or cuspy was still a matter of controversy till recent 
times. On one hand, many studies concordantly pointed 
to cored profile as a better fit to the rotation curve data 
("Burkcrt 19 951: iKravtsov et al. 1 11998*: ' Borriello fc Salucci I 
[2001 ; March csini et al. I l2002n Simon ct al. I l2005h . On the 
contrary, some further works seemed to contradict this 
conclusion with the fit to high resolution data being 
equ ally good for both cored and steeper profiles (se, 
e.g.. lSwaters et al. 112003 ). Although observational problems 
(such as slit offsets and beam smearing) and theoretical 
uncertainties (related to non circular motions and triax- 
iality) were originally thought to bias the determination 
of the inner slopes o f the density profiles from rotation 



curves d ata (see, e.g., Spekkens et al. 



ies (e.g.. lKuzio de Narav fc Kaufmann 



20051 ). recent stud- 
I2OIII ) have convinc- 



ingly shown that this is not the case so that one can safely 
rely on measured circular velocity profiles to investigate the 
cusp/core controversy. 

It is worth noting that cuspy haloes are predicted 
from dark matter only s imulations ([Navarro et al. 1996bl. 
I997I : iMoore et aLlliggsl : iPower et al. II2003I : iNavarro et al. I 



20041 . I2OI0I ) so that the impact of baryons, which actu- 



ally dominate the inner regions, is incorrectly neglected. 
It is interesting to note that the first solution proposed 
to the cusp /co r e probl em is connected to baryon physics. 
iNavarro et al. I (|l996al ) indeed studied how supernova - 
driven winds that expel a large fraction of baryons in 
the halo could give rise to a dark matter core. In the 
same framework, iGelman fc Sommer - Larsen I (|l999l ) repro- 
duced the observed rotation curve of DDO 154 by simu- 
lating NFW haloes and expo si ng th em to violent gas out- 
flows, while iGnedin fc Zhao I (|2002l ) considered the max- 
imal limit of complete blow out of all of the baryons. 
iRead fc Gilmore I (|2005J) then refined the analysis showing 
that, in order the previous mechanism to be efficient, it is 
necessary to have two impulsive mass -loss phases punctu- 
ated by gas re -accretion. More r ecent cosmological sim - 
ulations (JMashchenko et al. I I2OO6I : iGovernato et al. I I2OI0I ) 



^ For precision's sake, another type of LSB galaxies often dis- 
cussed in the literature are the massive, early -type, bulge - 
dominated LSB galaxies. These galaxies have properties entirely 
different from the late -ty pe LSB galaxies jSpravberrv et aL I 
ll995l:|Pickerine et al.]|l997l) . 



have continued to study the effec ts of baryons on t he in - 
ner DM inner profile. In particular, IGovernato et al. I (|2010l ) 
have recently showed that strong outflows from supernovae 
explosions in dwarfs -like systems may remove the low an- 
gular momentum gas inhibiting the formation of bulges 
and leading to shallow central dark matter profiles. Such 
a mechanism should, however, be inefficient in systems with 
a deeper potential well such as LSBs where the cusp/core 
problem is nevertheless present with several wo rks arguing 
in favour of both cored (|de Blok fc Bosma II2002I ) and cuspy 
(|van den Bosch et al. II200CI ) profiles. This circumstance has 
prompted the interest for alternative solutions to rescue the 
ACDM paradigm without drastic changes to its physics. 
An incompl ete list contains interactions of the DM with a 
stella r bar (|Weinberg fc Katz 1 12002| : iMcMillan fc Dehnen I 
l2005h. decay of binary black hol e orbits after galaxies merge 
( Milosavlievic fc Merritt I l200l|). baryon energy feedback 
from active galactic nucleus (|Peirani et al. 1120081 ). dynam- 
ical f riction of stellar/DM clumps against the background 
halo JEI Zant et al. ll2o"oIl . |2004 iRomano-Diaz et al. llJOOJI . 
l2009t ). 

A possible way out of the cusp/core problem is hidden 
in the recent result of de Blok et al. (2008) . Fitting the high 
reso lution rotation curv e data of the THINGS galaxy sam- 
ple jWaher et al. II2008I ). Ide Blok et ahl l|2008h have found 
that a NFW profile or a cored isothermal model fit equally 
well the data for galaxies with AIb ^ —19, while cuspy 
profiles are strongly rejected for systems with A'Ib ^ —19. 
In other terms, for low mass galaxies, a core dominated 
halo is clearly preferred over a cusp -like one, while both 
models are equally allowed for massive disk dominated sys- 
tems. That the density profile of dark matter haloes is 
not universal has been actually also sugg e sted by recent 



20001: iRicotti I l2003l: 
' 2009 : iMerritt et aTl 



numerical sinmlations (IJing fc Suto 

Navarro ct al. 2004 : [ Graham et al . . , . 

20061 : INavarro et al. l2010h showing that departures from the 
cano nical NFW p r ofile at 1 - 1 0% of the virial radius. More- 
over, iDel Popolo I (|2009l , |201G| ) used an analytical approach 
to show that such departures are stronger when baryons are 
taken into account with the inner slopes of the density profile 
becoming mass and redshift dependent (lower mass haloes 
having shallower profiles, with the slope becoming steeper 
towards higher z). 

Abandoning the universality of the halo density profile 
opens up the hunt for what is responsible for the variety of 
the halo properties which are inferred from observations. In 
order to shed light on this issue, it is worth looking at the 
statistical distribution of halo parameters and dark matter 
related quantities also looking for correlations with both the 
stellar properties and the environment the halo is embedded 
in. Motivated by this consideration, we have here analyzed 
the rotation curves of a sample of 27 dwarfs and 10 LSB 
galaxies collected from the literature. All of them have well 
measured rotation curves with both HI and H« data extend- 
ing to large R/Rd values (with Rd the disk scalelength) so 
that one is able to probe the dark matter dominated regions. 
Moreover, the stellar contribution to the circular velocity 
is well constrained by the measured luminosity profile and 
the constraints on the stellar population making it possible 
to constrain the stellar mass - to - light ratio and hence the 
baryon matter distribution. 

The scheme of the paper is as follows. The halo model 
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parametrization and the method used to constraint its pa- 
rameters are presented in Sect. 2, while the results of the ro- 
tation curve fits are discussed in Sect. 3. We then investigate 
the correlation of different halo properties with environment 
indicators in Sect. 4 and devote Sect. 5 to conclusions. 



2 DATA AND METHOD 

Being, by definition, not observable, the dark matter haloes 
can be probed only indirectly by their effect on the galaxy 
dynamics or through their lensing effect. Fitting the rota- 
tion curve is the only way to constrain the properties of 
the dark matter component for local dwarfs. To this end, 
we have therefore searched the literature for systems with 
high quality rotation curve data probing the gravitational 
potential with good sampling, large radial extension (i.e., 
up to i? > Ropt = 3.2Rd with Rd the disk scale length) and 
small errors. Whenever possible, we rely on Hq data or a 
combination of HI and Ha so that the potential impact of 
beam smearing on the circular velocity is reduced and does 
not bias the estimate of the inner slopqj. The data typically 
provide the total circular velocity so that one needs first 
to model the baryons (stars and gas) contribution. While 
this is not a problem for the gas (whose mass profile is di- 
rectly obtained from HI data), an estimate of the stellar 
mass - to - light ratio T* is needed in order to convert the 
observed luminosity profile into a mass one. Should galaxy 
colors be available, one could rely on stellar population syn- 
thesis models to infer T^ which is indeed the case for most 
of the galaxies in our sample. On the contrary, when this 
is not the case, a fiducial value can be adopted based on 
the prior constraints on T*. We do not revisit here this 
issue setting T* to the value adopted in the paper which 
the data are presented in. However, in order to take care of 
possible deviations from this value (due to possibly wrong 
assumptions or uncertainties in the different quantities en- 
tering stellar population synthesis codes), we add a free pa- 
rameter «:* = Tt/T fid with T* the actual unknown stel- 
lar M/L ratio and T fid the assumed one in the original 
paper. Note that k* can also be smaller than unity, a typ- 
ical example being the assumption of an incorrect initial 
mass function (scaling down the M/L estimates by a con- 
stant factor smaller than unity). We therefore conservatively 
set a fiat prior on k, over the range (0.55, f .45) so that we 
fully take into account both the systematic uncertainties on 
the initial IMF choice and the scatter introduced by vari- 
ations in metallicity and star formation rate. Although we 
are mainly interested in dwarfs, we initially select also LSB 
galaxies since they share most properties with the dwarfs 
ones. We finally end up with a sample of 41 galaxies reduc- 
ing to 37 (28 dwarfs and 9 LSBs) after removing multiple 
cases. It is, indeed, sometime possible that the same galaxy 



^ A possible concern on the use of the Ha data is related to H« 
possibly tracing biased star forming regions of the disc (as, e.g., 
spiral arms). As a consequence, the measured circular velocity 
could be non representative of the disc stars kinematics. While 
we can not definitively exclude such a possibility, it is nevertheless 
worth noting that, whenever both HI and Ho data are available 
in the same region, they closely track each other. We are therefore 
confident that no significant bias is induced by the use of Ha data 



has multiple measurement of the rotation curve in which 
case we choose the better quality one. We have, however, 
checked that the results fitting the other dataset are consis- 
tent thus making us confident that no bias has been induced 
by this choice. Table 1 gives the list of galaxies and the ref- 
erences for the data on the rotation curve. In more detail, 
dwa rfs galaxies main ly come from the lSwaters et al. I (|201ll ) 
and lOh et al. I (|201l|) samples and have a typical dynamical 
mass ~ 10^ M0. Another dwarf spiral (NGC 2976) comes 
from the THINGS sample which also gives us two dwarfs 
irregular (IC 2574 and NGC 2366). Finally, the sample in 
ISimon et all (|2005l ') contributes the dwarf spiral NGC 4605, 
the SBc galaxy NGC 5649 and the Sc one NGC 5963, and 
the low mass spiral NGC 6689 with different listings for its 
morphological type. LSB s systems mainly comes from the 
Ide Blok fc Bosma I (|2002l ) sample which we refer the reader 
for details. 

In order to fit these data, we have first to choose a den- 
sity profile for the dark matter distribution. To this end, one 
can rely on the outcome of N - body simulations choosing a 
double power-law model, i.e. p oc x~°'^"{l + x)~"°"'^^"'"' , 
with X = r/Rs and (oin, ciout, Rs) setting the shape of the 
profile. While there is a large consensus on aout = 3, a still 
open debate exists on which ai„ value better fits the nu- 
merical simulation s and whether this i s or not a universa l 
quantity (see, e.g.. Navarro et al. 1 1997 : Moore et al. 1119981 : 



iJing fc Suto l2000l : lRicotti l2003l : lDiemand et al. |2004| ). Ac- 
tually, more recent simulations suggests that the logarithmic 
density slope a = dlnpoAi/dlnr never asymptotes a con- 
stant value in the small r limit, while an e xponential - like 
decline is achieved in the outskirts. FoUowina lNavarro et al. I 
(1200J), we therefore assume that the halo density profile is 
given by the Einasto (1965, 1969, see also [ Cardone et al. I 
I2OO5I : iMamon fc Lokas II2005I : iGraham et al. II2OO6I ') model: 



PDAiir) = p-2 exp < — 2ni3j\ 



R- 



(1) 



with udm a slope parameter, R-2 a characteristic radius 
defined so that a{R-2) = —2 and p-2 a scaling den- 
sity. It is convenient to reparameterize the model replac- 
ing the {R-2, P-2) parameters with {c^ir,M^ir) where we 
have introduced the concentration Cvtr ~ R-2/Rvir, the 
virial mass Mvir (defined as the mass within a sphere con- 
taining A.uir — 337 times the cosmic mean matter den- 
sity pM = 3HqQ,m/8itG) and the virial radius R^ir — 
{SMvir/^TvAvirPMy. Actually, rather than Mvir, we will 
use as a parameter 14ir = {GM-uir / RvirY , i-e. the circular 
velocity at the virial radius. Assuming spherical symmetry, 
the mass profile simply reads : 



M{r) = 



^TmDMR-2P'2e 

(2noA/)^"oA' 



!2"DM 



X T|3no„,-^ ) (2) 



with 'y{3nr)M,x) the incomplete gamma function. The cir- 
cular velocity is then v^{r) = GM{r)/r and may be analyt- 
ically estimated as function of the the model parameters. 

As already hinted at above, we parameterize the Einasto 
model through the three quantities {nDM,Cvir,Vvir) to 
which we add k* to take care of the uncertainties on the 
stellar M/L ratio. In order to determine the four quantities 
(k*, TiDA/, c„i,., Kiir), we first define the likelihood function 
£(p) oc exp hx'(p)/2] with 
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Table 1. Galaxy sample and constraints on model parameters. Columns are as follows: 1. galaxy id; 2. type (0 for dwarfs, 1 for LSBs); 
3. reference for the rotation curve data according to the following scheme: Oh = Oh ct al. (2011), dBB02 = de Blok & Bosnia (1992), 
THINGS = de Blok et al. (2008), SOS = Simon et al. (2005), Swll = Swatcrs et al. (2011); 4. tidal index; 5. projected distance of the 
tenth closest galaxy (in kpc); 6. number of galaxies within a projected distance of 750 kpc; 7. - 11. median and 68% confidence range for 
the model parameters (T/Tj-j^,a, c„iri l^uin logM^i^) assuming an Einasto density profile for the DM halo; 12. reduced x^ = x^ jd.o.f. 
for the best fit model. Note that, for galaxies in the S05 sample, we have to set the stellar M/L ratio to its fiducial value since we have 
only data on the DM circular velocity so that no constraints on T*/Tj^irf = 1.0 can be reported in the table. 



Id 


Type 


Ref 


e 


Rw 


A/Vso 


T*/Ty,d 


riDM 


^vir 


K>r 


log Myir 


-2 

X 


DD053 





Oh 


0.7 


411.724 


24 


0.811-- 


129+^-" 
-"-■^■^-0.83 


s-^tii 


27.5lf3:« 


W.otol 


0.19 


DD0185 





dBB02 




342.667 


26 


0.80+i;i^3 


63+^" 

'-'■"•^-0.46 


ii-5+^:i 


56.11^8^ 


io.9+i-;g 


1.45 


Hoi 





Oh 


1.5 


187.429 


27 


U./0_0.i4 


n 90+009 

'-'■^^-0.08 


18.6+°-? 


4+0.2 
8.4_o.2 


S 4+0-1 
8.4_o.l 


2.18 


HoII 





Oh 


0.6 


445.606 


17 


n 74+0-29 


12 4+4-3 


3.6+2-« 


29.8lt.} 


10.lt°-2 


2.37 


IC2233 


1 


dBB02 




20.943 


8 


Qg-^ + 0.37 


2-00+J-^^ 


7-i^^:^ 


134.2l^«-* 


12 l+i-'i 


0.70 


IC2574 





THINGS 


0.9 


172.046 


26 


7fi+0-20 
'J- '0-0.14 


-, 17-1-0.73 


6.3t?i 


62.9t^I:i 


n i+°-8 


0.50 


M81dwB 





Oh 


-0.9 


119.732 


21 


0.86+«-g 


2 64+^-** 


24 9+^-9 

^"'•y-14.3 


17.4l?2^5 


04+1.4 


0.08 


NGC2366 





THINGS 


1.0 


525.228 


10 


0.84+«;g 


1 04+°-*^ 
^■^*-0.30 


12.8+;-^ 


21.9+^-^ 


Q 7+0.3 

^.'-0.2 


0.50 


NGC2976 





THINGS 


2.7 


131.728 


27 


0.701°:?? 


O-SOtgaf 


io.8i°;« 


69.0lf3'o' 


ii.2tJ:^9 


0.50 


NGC3274 


1 


dBB02 


-0.3 


310.581 


14 


n 77-1-0.38 
^•' '-0.17 


9 nH+O-57 

2-08-0.42 


28.6+2-1 


29.4t^i 


io.it°-^ 


0.49 


NGC4605 





S05 


-1.1 


417.435 


18 




60+°" 


28.5t°;^ 


19.9+11 


9.6t;j:l 


1.37 


NGC5023 


1 


dBB02 


-0.5 


210.311 


24 


76+°-8'' 


1 oo+°-''2 

^■^'J-0.30 


18.3+^° 


29.5ttl 


io.it°:? 


0.19 


NGC5949 





S05 




479.538 


3 




3-98t^:f8 


8.9t^-.°3 


129.9+^^°^^ 


i2.iti:^ 


0.25 


NGC5963 





S05 






314.317 


6 




9 1+^-2 


4 7+8.9 

^•'-3.3 


222.6+12V 


12 7+°-5 


0.42 


NGC6689 





S05 








2 




3.27+2-22 


g 4+4-1 


150.6+278^ 


12 2+1-2 


0.31 


UGC731 





Swll 








1 


O-SS+o'al 


3.79+2-99 


i5.8i?:« 


46.2is;2 


10 7+0.6 


0.11 


UGC1230 


1 


dBB02 






142.162 


2 


0-83+°-^I 


i-8o+°:«5j 


14.31- 


52.7l«-? 


10.8+°-? 


0.50 


UGC1281 


1 


dBB02 


-1 


.2 


324.170 


17 


PI 7-1+0:21 


0.54tg-jg 


ii.4i;-f 


22.6IJI52 


9.71°:^ 


0.07 


UGC3137 


1 


dBB02 






429.129 





61+°" 


^9+0-10 


9+0.2 

8-^-0.3 


49.3t°:^ 


io.7tg:} 


2.92 


UGC3371 


1 


dBB02 






517.837 


1 


^•'^-0.17 


1 40+0-86 


9 1+1.4 
y.±_-[ 7 


52.9tlli 


10.8+°-« 


0.01 


UGC4173 





dBB02 






587.408 


1 


84+0-3'^ 


2 0-1+1.42 


40+1.5 


42 9+34-3 


io.7t°;? 


0.06 


UGC4325 


1 


dBB02 






315.034 


4 


n S4+0-39 


i-i9tg:S 


16.0+^-9 


Sfi c+107.4 
86-8_36.7 


11 c + l.O 
ii.b_l37 


0.01 


UGC4499 





Swll 






374.473 


6 


0.83t°:^5 


2.58t?-.^° 


11.3l^;« 


44 2+3^-1 


10.6t°1 


0.33 


UGC5414 





Swll 






402.487 


1 


0.86l°:i^ 


2 00-1-1.93 


9-^tll 


50.3+21? 


io.8t°;? 


0.12 


UGC6446 





Swll 






173.865 


15 


0.96l°i^ 


90+2-40 
■^■^'J-0.93 


15.6+li 


42.9+^V 


10.6«-3 


0.22 


UGC7323 





Swll 






93.111 


67 


0.911^:^7 


9 -50+1-87 

^.OU_Q gg 


9.2tli 


76.7+Ji's' 


ii.stJ:? 


0.22 


UGC7399 





Swll 






87.709 


53 


1-031^:?^ 


4 84+1-3S 


20.8+^-^ 


65.4+i°i« 


ii.it^:^ 


1.55 


UGC7524 





Swll 





1 


160.137 


142 


n 00+0-36 
U-y^_0.28 


2 04+1-23 


14.0+°:^ 


40.0ti^g° 


10.5t°-* 


0.18 


UGC7559 





Swll 





1 


169.059 


83 


0.84t°:g 


2 9c:+3.01 


14.6+3:4 


19.lt?'4* 


q c + 1.2 


0.03 


UGC7603 





Swll 






228.234 


144 


o.95t°:^g 


-, 09+1.13 
-■-■8^-0.66 


14.61^;^ 


29.9+^7^ 


IO.I+0I 


0.17 


UGC8490 





Swll 


-1 


.1 


413.572 


10 


0.83t°:i? 


4.08tJ;2^ 


25.512:2 


io.rtll 


io.5t°-:? 


0.11 


UGC9211 





Swll 






395.074 


2 


n 00+0.34 


2 29+3.45 
^■^^-0.99 


14-41^:^ 


32.9tl%° 


10.2t°:«3 


0.03 


UGC10310 


1 


dBB02 











0.82l°:33 


1 71+l-''2 
l-'^-0.93 


9+2.1 


4Q 7+61-1 


io.8ti:°7 


0.09 


UGC11707 





Swll 








3 


0.83+°-^? 


c CK+4.74 


9+4.0 
8.2_4.5 


09 c+20.4 
82.O-24.0 
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x^(p) = E 



VcobsRi - Vc,th{Ri,p) 
Si 



where the sum is over the A/'ots datapoints with Vc,obs{Ri) 
and Ei the i- th measured value and its error and Wc,th(-Ri, p) 
the theoreticaUy expected one. To efficiently explore the pa- 
rameter space, we use a Markov Chain Monte Carlo method 
running three chains and checking their convergenc e through 
the Gellman-Rubin test lCelman fc Rubin I (|l992l ). Accord- 
ing to the Bayesian statistical theory, the set of parameters 
maximizing the likelihood do not necessarily represent the 
most confident constraints on each single parameter. On the 
contrary, one has to look at the marginalized distributions 



and take the median and 68% confidence ranges as final con- 
straints. To this end, we first cut out the burn -in phase of 
each single chain, then merge them together and use the 
resulting sample (after thinning to avoid spurious correla- 
tions) to estimate the median and 68% confidence ranges. 
This sample can also be used to infer constraints on derived 
quantities (such as the virial mass and the logarithmic slope 
Q at different radii) by estimating them along the merged 
chain and studying the distribution of the resulting values. 
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Figure 1. Best fit curves superimposed to tiie data for tiiree typical cases with very low (left), average (centre) and high (right) x^ . 
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Figure 2. Scatter plot of the MCMC samples for the fit to the UGC 11707 data. 



3 RESULTS 

As a general remark, we find that the Einasto model in- 
deed fits quite well the rotation curve data for both dwarfs 
and LSB galaxies so that the estimated constraints on the 
model parameters (reported in Table 1) may be deemed as 
reliablqj. Table[T] also gives the values of the reduced x^ , 
i.e. x^ = X^ /d-o.f. with d.o.f. — Afots — 4 the number of 
degrees of freedom. Actually, this quantity should be con- 
sidered only as a way to quantitatively guess the accordance 
of the model with the data, but one must not interpret it 
as the textbook x^ estimator hence expecting X^ ^ 1 for 
the best fit model. Indeed, the errors on the measured cir- 
cular velocities are not Gaussian distributed and, moreover, 
typically takes into account also corrections for non circu- 
lar motions and asymmetry between the approaching and 
receding side of the rotation curve data. Finally, for some 
cases (as, e.g., the de Blok & Bosma (2002) sample), the 
data have been suitably smoothed so that a good model 
will fit the data with a very small x^- This is indeed the 
case for many curves as can be seen, e.g., for UGC 4325 
plotted in the left panel of Fig.[T] While this is an overopti- 
mistic case, the agreement within the theoretical prediction 
and observed rotation curve is typically quite good with the 
best fit curve closely interpolated the measured circular ve- 
locity (as, e.g., for UGC 11707 shown in the central panel 
of Fig.[l]). For few cases, we get a large x^ which could be 
interpreted as a failure of the model. Actually, rather than 
signalling a mismatch between data and theory, a large x^ is 
typically an evidence for the presence of wiggles in the data 



^ In the appendix, we will briefly present an analysis of simu- 
lated rotation curves to further check that the fitting procedure 
is indeed able to correctly recover the input model parameters. 



due to residual non circular motions or clumpy gas distribu- 
tion as can be seen, for instance, in the case of Holmberg I 
in the right panel of Fig.[T] We therefore conclude that the 
Einasto model provide a good description of the DM haloes 
of both dwarfs and LSB galaxies and rely on it to estimate 
the quantities of interest for the following analysis. 

As it is apparent from the 68% confidence ranges re- 
ported in Table[T] the constraints on the fitting parameters 
are not quite stringent. Needless to say, this is a consequence 
of both the data quality and the degeneracy in the 4D pa- 
rameter space. Indeed, in order to break the degeneracies, 
one should not only trace the circular velocity up to large 
R/Rd values (as is the case for most galaxies in our sample), 
but also reduce as much as possible the measurement uncer- 
tainties and finely sample the rotation curve. Indeed, the 
smaller errors are ob tained for the galaxies in the THINGS 
(|de Blok et al. Il200a ) sample which are close to fulfill these 
requirements. Actually, also in these cases, the constraints 
on the virial velocity and hence the virial mass can be weak 
if the data do not probe the outer halo dominated regions. 
As a consequence, the stellar M/L ratio becomes degen- 
erate with log Mvir since both quantities scale the relative 
contribution of the stars and DM to the total Vc{r). Lack- 
ing a strong determination of T*, one can not weight the 
DM contribution in the disk dominated regions so that the 
uncertainties on the slope noAi and concentration Cvir in- 
crease. Notwithstanding these problems, it is nevertheless 
interesting to look at the distribution of the halo parame- 
ters {udm, Ciiir, log Mvir) where here and in the following we 
have replaced the virial velocity with the logarithm of the 
virial mass since this latter quantity provides a better char- 
acterization of the DM halo. Moreover, we impose a loose cut 
on (tida/, Ci.ir, log Afuir) excluding galaxies with udm > 11 
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or log M„ir < 9.0 since they are likely outliers thus ending 
up with a sample made by 27 dwarfs and 8 LSBs. 

As a first test, we compare the distribution of the model 
parameters for dwarfs and LSBs. The mean, median and 
68% confidence ranges of (udm, c^tr, log Mvir) turn to bqj: 

{riDM) = 3.05 , fiDM = 2.58 , 68% CL = (1.29, 4.36) , 

{c^ir) = 13.51 , c^ir = 12.76 , 68% CL = (8.14,20.81) , 

{fivir) = 10.71 , jlvir = 10.59 , 68% CL = (9.99, 11.42) , 

for the dwarfs sample (with fj,vir = logM^ir) and : 

{riDM) = 1.48 , fiDM = 1.49 , 68% CL = (1.00, 2.00) , 

(c„i^) = 14.13 , c^ir = 11.39 , 68% CL = (8.23, 18.28) , 

{fivir) = 10.73 , fivir = 10.76 , 68% CL = (10.07, 11.48) , 

for the LSBs sample. Although the halo mass range probed 
is almost the same, the slope parameter udai takes larger 
values for the dwarfs than for the LSBs. This is essentially 
due to the strongly asymmetric udai distribution for the 
dwarfs having a long tail towards high udm- As a conse- 
quence, both the mean and the median are larger than for 
the LSBs. On the contrary, we find no significative difference 
for the concentration parameter with the 68% confidence 
range being almost perfectly overlapped. Actually, the sta- 
tistical significance of the udm difference is undermined by 
the large errors on the individual determinations. Moreover, 
the statistical properties of the halo parameters distribution 
for the joint sample, namely 

{udm) = 2.69 , fiDM = 2.22 , 68% CL = (1.17, 4.08) , 

(c„i,.) = 13.65 , c^ir = 12.76 , 68% CL = (8.14,20.81) , 

{^JL^ir) = 10.71 , fi^ir = 10.60 , 68% CL = (9.99, 11.48) , 

are quite similar to the dwarfs one so that we prefer to in- 
crease the statistics adding LSBs to the dwarfs sample. Al- 
though such a strategy decreases the homogeneity of the 
sample, we have checked that it does not bias the results of 
interest as could be anticipated noting that the final sample 
is still dominated by dwarfs. In the following, for sake of 
simplicity, we will refer to the full set as the dwarfs sample. 
Since N - body simulations point at the Einasto model 
as the best representation of the dark haloes density pro- 
file, the successful fits of the dwarfs galaxies rotation curves 
may be seen as a further evidence of this model virtues. 
Actually, our result tells us something more since it refers 
to a mass range typically not resolved by simulations. It is 
therefore interesting to compare udm distribution we find 
with the ones predicted by the numerical simulations. These 
works find values for n pM covering th e range (5.0, 7.0) in a 
roughly uniform way ([Navarro et al. 1 2004: Merritt et al. I 



I2OO6I : iGraham et al. II2OO6I : iNavarro et al. Il2010h . while our 
Udm distribution, although overlapping with the predicted 
one, is definitely more populated by shallower profiles. This 
can also be read from Fig. [3] where we plot udm vs logM^ir 
for both our data and the results from N - b ody simulations 
iNavarro et al. 1 12004| : iMerritt et al. I I2OO6I '). A straightfor- 
ward comparison is, however, not possible given both the 
different mass range and the uncertain impact of baryons 

^ Hereafter, we denote with [x) (£) the mean (median) x value. 




Figure 3. Einasto slope parameter nuM vs the virial mass 
log M^ir for our sample (black) and the N - body haloes of Navarro 
et al. (2004, blue) and Merritt et al. (2006, red). 



(which are not accounted for in the simulations and may 
change the density slope). In order to partially alleviate 
this problem, it is much more instructive to compare our 
noM va lues to those o b taine d by fitting the full THINGS 
sample. IChemin et al. I l|201ll 'l have indeed found that the 
Einasto model fits the observed rotation curves in a very 
good way being also superior to popular models as the NEW 
and the isothermal ones. Moreover, they find that better 
fits are obtained by adopting a Kroupa (2001) IME rather 
than the fiducial diet Salpeter used by the THINGS team. 
This is the same as downscaling the stellar M/L by a factor 
~ 1/1.16 ~ 0.86 which is co nsistent with the typ ical values 
we find for our k* parameter. IChemin et al. I ()201lh have also 
found that their halo sample can be divided in four classes 
according to the udm value with the most numerous class 
being represented by systems with 0.1 ^ tidm ^ 4.0 and 
the rest of galaxies mainly having udm > 4.0 values. This 
is roughly consistent with what we find here thus suggesting 
that dwarfs and normal spiral galaxies cover the same range 
for the slope of the Einasto model. However, since our con- 
fidence ranges are quite large and strong degeneracies with 
the other model parameters are present (as can be seen from 
the distribution of the MCMC points in Eig.[2]), we prefer to 
not speculate further on this issue, but only point out at this 
qualitative agreement as a further evidence of the reliability 
of our choice of the Einasto model for the dwarfs galaxies 
haloes. 



4 DARK HALO PROPERTIES 

Having convincingly demonstrated that the Einasto model 
provides a good fit to the rotation curve data, we can use 
the constraints derived on its parameters to estimate other 
quantities of interest. To this end, one has only to evaluate 
y = yijiDM, Ctiir, log Mvir) aloug the final merged chain and 
then look at the distribution of the y values. It is then in- 
structive to look for correlations of the estimated quantities 
with both the stellar luminosity and size and the environ- 
ment the dwarf resides in. To this regard, it is worth stressing 
that our aim here is not to constrain any particular galaxy 
formation scenario, but mainly to search for any evidence of 
scaling relation which could then be used as a guidance to 
improve our understanding of haloes assembly. 
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Table 2. Correlation properties of the Einasto model parameters with the environment. Rows are as follows: 1. Spearman correlation 
coefficient; 2., 3., 4. median and 68% CL of (a, b, cri„t)- The number of galaxies used is A^ = 13, 29, 33 for correlations with 0, Eio, Dtsq, 
respectively. See text for details on the fitted quantities and the median values of the input parameters. 
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4.1 Preliminary remarks 

Before discussing the results of our analysis, we have first 
to explain what are the estimator we use to describe the 
environment a galaxy lives in. To this end, we use three 
different statistics to quantify whether a galaxy lives in a 
dense environment or is rather an isolated one. First, we 
consider the tidal index (JKarachentsev fc Makarov Ill996l ) 



e = max[log(Mfc/DE]+C 



(3) 



where {Mk,Dk) are the mass and distance to the chosen 
galaxy of the fc-th system and C is a normalization con- 
stant. Estimating Q for a given galaxy is actually quite dif- 
ficult and typically possible only for systems in the Local 
Group. A cross match be tween our sample and the Catalog 
of Neighboring Galaxies (jKarachentsev et al. 1120041 ) allows 
us to infer for 10 dwarfs and 4 LSBs only. 

In order to further test a possible correlation between 
the inner slope and the environment, we resort to two other 
commonly used indicators. The first, denoted as Eio, is the 
projected number density of galaxies inside the circle of 
radius equal to the projecte d distance of the tenth near- 
est neighbor ({Dressier lllQSOf ). while the second, denoted as 
D750, is the projected number density of galaxies within a 
circle of radius equal to 7 50 kpc centred on the system of 
interest (JTreu et al. II2009I ). We will refer to Eio and D750 as 
the local and global overdensity and stress that both these 
quantities are normalized with respect to their mean value 
over the sample so that galaxies with Eio and -D750 larger 
(smaller) than unity are in over (under) dense regions. As 
a general remark, we note that, although easy to calculate, 
both Sio and D750 are looser indicators of dynamical state 
of a galaxy with respect to the tidal index. For instance, a 
galaxy could have a large value of Eio, but a small O if it is 
embedded in an environment rich of small mass systems. In 
such a case, the galaxy should be considered isolated even 
if Eio is large. However, such case should be rare and, in- 
deed, we find that both Eio and D750 positively correlate 
with O so that we use them as supplementary tools to check 
whether the environment impacts the halo density profile. 
Using the NED databasqj, we are able to estimate Eio for 
30 galaxies (22 dwarfs and 8 LSB), while -D750 is available 
for 34 systems (26 dwarfs and 8 LSBs) thus doubling the 
sample usable for the correlations with the tidal index. 
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To quantify whether there exists any correlation be- 
tween a given quantity yinoM, c^ir, log M^tr) correlates and 
the environment indicators (O, Eio, D750), we use two dif- 
ferent approaches. First, in order to get a quick result, we 
compute the Spearman rank correlation coefficient denoted 
as C{x,y). Since the sample used is not large enough, it is 
possible that few points drive the C{x, y) estimate. In order 
to avoid such a bias, we resort to a jacknife approach eval- 
uating C{x, y) for A/" — 1 samples obtained by excluding one 
galaxy at time and use the mean value and the standard 
deviation of this sample as our final estimate. To quantify 
whether the result indicates the presence or lack of a statis- 
tically significant correlation, we can compute the z- score 
defined as 2 = [(A/" - 3)/1.06]^''^arctanh[C(x,y)] where A/" is 
the number of datapoints. Since z approximately follows a 
standard normal distribution, one can quickly check whether 
the C{x,y) value corresponds to a statistically meaningful 
correlation. 

The Spearman rank estimator does not take into ac- 
count the errors on y coming out from the uncertainties on 
the constrained halo parameters. We therefore rely on a dif- 
ferent approach briefiy sketched below. Let us suppose that a 
linear relation exists among two different quantities denoted 
(R, Q) such as R = aQ + b and denote with aint the intrin- 
sic scatter around this relation. Calibrating such a relation 
means determining the two coefficients (a, b) and the intrin- 
sic scatter aint- T o this aim, we will re sort to a Bayesian mo- 
tivated technique ID' Agostini I (|2005l ) thus maximizing the 
likelihood function £(a, 6, atnt) = exp [—L{a, b, aint)] with ; 



L{a,b,aint) 



+ 



1 -^ _{R, - aQ, - &)2 
2^ a?. 



j2 

'mt 



(4) 



where the sum is over the J\f objects in the sample. Note 
that, actually, this maximization is performed in the two 
parameter space {a, aint) since b may be estimated analyti- 
cally as : 



Ri — aQi 






<^int + ^R, + (^ <^t 






'Q. 



SO that we will not consider it anymore as a fit parameter. 
Two caveats are in order here. First, as we will see later, 
we are actually looking for power -law like relations so that, 
when needed, we will use logarithmic units in order to lin- 
earize them. As a second issue, we note that Eq.Q implic- 
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Figure 4. Median normalized Einasto slope parameter nuA/ vs the three environment indicators. 
Table 3. Same as Table[2]but for the logarithmic slopes (a^, qq.i) using Af = 12, 25, 28 galaxies for correlations with (©, Sio, -D750). 
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itly assumes that the error on the quantities involved are 
equals on both the positive and negative side, i.e., the 68% 
confidence range is symmetric around the mean. This is not 
the case for m ost of the quantitie s we are interested in so 
that, following Id ' Agostini I (|200J), we take the mean of the 
positive and negative errors as input to the likelihood func- 
tion. Moreover, we typically fit log {y/y) rather than logy in 
order to reduce the degeneracy among the slope a and in- 
trinsic scatter aint- We finally look at the median and 68% 
confidence range of the slope a to check whether the two 
quantities are correlated or not noting that a null value is a 
definitive evidence for a lack of any correlation. 



4.2 Halo model parameters 

We start investigating correlations between the environment 
and the halo model parameters {nr>M,Cvir, fJ-vir) where, 
hereafter, we will use the shortened notation ^^ir to refer to 
log M^ir- The results obtained are summarized in Table[2l 

Let us first consider the correlations involving the 
Einasto slope parameter udm which we fit as a functioiij 
of (9,logEio,logD75o) using as variable log{nDM/nDM) 
with hoM = 2.22 the median value. Both the Spearman 
correlation coefficient C and the slope a of the fitted rela- 
tions strongly indicate that udm does not correlate with the 
environment. Actually, although we refer to it as the slope 
parameter of the Einasto profile, the logarithmic density is 
not related to udm only, but rather to the full set of halo 
parameters being 



^ Here and in the rest of the paper, we will use logarithmic vari- 
ables for (Sio,i?75o) so that correlations with this parameter 
will be of the form y ex x°' with y the quantity of interest and 
X = (Sio, Drso)- On the contrary, we will retain linear units for 
© since the tidal index is already a logarithmic quantity. 



«(0 - -2 ( ^ 



= -2 






(5) 



The slope of the density profile better characterized by aivs) 
where Ta is some reference radius. As first possible choice, 
we set Ts = Rd since the region around the disk scalelength 
is typically the best sampled by the rotation curve data. 
However, since Rd is different from one case to another, we 
also use rs = O.lRvir so that we sample a at the same 
scaled radius for all the galaxies. Moreover, as can be easily 
understood from Eq.©, the smaller is rs/Rvir, the better 
is oi{rs) constrained so that setting rs/Rvir = 0.1 represents 
a good compromise between reducing the uncertainties and 
avoiding to infer constraints on a quantity defined over a 
range not probed by the data. In the following, we will refer 
to a{Rd) and a{0.1Rvir) as ad and ao.i, respectively. Note 
that Od is evaluated at a smaller radius than ao.i given 
that, in order to have O.lRvir/Rd > 1, it is sufficient that 
Rvir/Rd > 10 holds which is always the case. 

The mean, median and 68% confidence range reaqj: 

(Qd) = -1.03 , &d = -1.02 , 68% CL = (-1.52, -0.51) , 

(qo.i) = -2.28 , ao.i = -2.33 , 68% CL = (-2.63, -1.76) , 

with both the distributions being almost symmetric around 
the mean values. Considering that |a(r)| approaches the null 
value in the very inner regions (r « R-2) and that, for a 
typical Rvir/Rd ~ 60 value, ao.i refers to the slope of the 
density profile in the terminal part of the rotation curve, it 
is easy to qualitatively understand why the Einasto model 
turns out to be so successful in fitting the data. Indeed, the 



^ Note that we exclude from this analysis the four galaxies from 
the Simon et al. (2005) sample since we do not have an estimate 



of their disc scalelength. 
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Figure 5. Same as Fig. [3] but for the DM mass fraction at the disc scalelength. 



Table 4. Same as Table[2]but for the DM mass content {fj, /-2) using J\f = 11, 25, 28 galaxies for correlations with (0, Eio, i?75o)- 
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dark haloes has a very shallow density profile in the disk 
dominated region (r << Rd) thus mimicking cored models 
which are well known to fit the dwarfs galaxies rotation curve 
data. In the intermediate r ~ Rd regions, the logar i thmic 
slope takes similar values to the one of the iBurkert I (|l995l ) 
model which is known to work well in this regime. Finally, 
the Einasto model mimics the isothermal one (which has a 
constant logarithmic slope a = —2) over the r » Rd range 
probed by the data. It is, therefore, this ability to interpolate 
among different models which makes the Einasto profile so 
successful in reproducing the rotation curve data of both 
dwarfs and normal spiral galaxies. 

The quite large confidence ranges for (n_Di\/, a^, ao.i) 
may also be considered as an evidence against the existence 
of a universal halo profile. It is then worth investigating 
whether the environment takes some role in driving this non 
universality. Should this be the case, we must find a corre- 
lation between the environment and slope indicators. The 
correlation coefficient and the {a,b,(Jint) parameters of the 
fitted linear relations (using ai/oi with i — d, 0.1 as vari- 
ables) are summarized in Table[3l Although C{x,y) can be 
sometimes large (as for the ao.i - © case) or the slope a being 
not consistent with the null value within the 68% CL (as for 
the Ud-Drso case) but with a large scatter, we can safely 
conclude that there is not any statistical evidence of a cor- 
relation of the logarithmic slopes with the environment. As 
an example, we plot udm vs the environment indicators in 
Fig.[4] giving a visual confirmation of the former conclusion. 

The same conclusion is strongly suggested by the results 
in Table[2] for the concentration parameter and its correla- 
tion with the environment obtained by fitting log {Cvir/cvir) 
vs {&, log Eio, log D750). It is particularly instructive to note 
that, although C{Q,Cvir) is quite large hence arguing in 
favour of a strong correlation, the slope a, on the contrary, 
definitely shows that such a correlation does not exist at all. 



Such an example, therefore, highlights the importance of 
fully taking into account the errors on the quantities before 
drawing any conclusion on a relation between two quantities. 
Somewhat surprisingly, we find a strong correlation 
among the virial mass and the local and global density es- 
timators (Eio,D75o), but not with the tidal index O. Con- 
trary to what happens for Cvir, the C(x,y) values are al- 
ways quite small, but the slopes a of the correlation^ with 
(logEio, logD75o) turn out to be non vanishing at the 68% 
CL thus arguing in favour of a dependence of /iui,- on the 
environment (although with a non negligible scatter). It is 
worth noting that both Eio and D750 are not related to 
the dynamical state of the environment and can be contam- 
inated by projection effects. As such, one can not exclude 
that the found correlation is only an artifact. Should we take 
them for actual relations, we should conclude that the higher 
is the local density Eio, the less massive is the dwarf halo. 
On the contrary, the higher is the global density D750, the 
larger is the halo mass. Reconciling these two somewhat op- 
posite behaviour is actually quite difficult so that we warn 
the reader to not overrate the significance of the detected 
^vir - Eio and ^vir - D750 correlations. 



4.3 Dark to stellar ratios 

In order to shed some light on the interplay between the 
baryonic component and the dark halo of galaxies, it is in- 
teresting to investigate how they contribute to the total mat- 
ter distribution. We have therefore estimated the DM mass 
fraction /njif (r) = MDM{r)/[Mbar{r) + MoAf (r)] where we 



define an effective total baryon mass as h'harir) 



rT/G 



* Since fimr is already a logarithmic variable, one could use 
l^vir / (J-uir as fitted quantity. We rather prefer to use fJ-vir — jj-vir = 
log {Myir / Mvir) with M^ir = dex(/i„ir) — 3.98 X 10*° Mq. 
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Table 5. Same as Table[2]but for the size ratios (7^—2, 'R-vir) using J\f = 11, 25, 28 galaxies for correlations with (©, Eio, O750). 



Id 


7^_2-e 


7^-2 -Sio 


■R-2-Dr5o 


IZyir - 


T^uir - Sio 


Tlvir - D7S0 


C{x,y) 


0.00 ±0.14 


-0.13 ±0.05 


-0.21 ±0.04 


-0.34 ±0.14 


-0.25 ±0.05 


-0.30 ±0.03 


a 


-o.oatH;; 


-0 06+°-^^ 

U.UD_Q 21 


n r,/|+0.15 


-0.031°;°^ 


_Q 14+0.28 
'^■-'-+-0.24 


'^■^'^-0.25 


b 


-0.09tH? 


—0 09+0.36 


0.00l°;f7 


-0.03t«:«l 


„ 2-1+0.35 


1Q+0-31 
'^■^^-0.29 


O'int 


omtZl 


n lfi+0-05 


n ifi+O'O'' 

^•^O„Q04 


0.24t-.°^ 


24+°-''5 


24+°-''5 



Table 6. Same as Tablo[2]but for the Newtonian accelerations at Rj using A^ = 12, 26, 29 galaxies for correlations with (0, Sio, Dyso)- 



Id 


^f«-e 


7r^-sio 


l^^'-Dno 


7r-e 


-Td^'-SlO 


tS^'-^^tso 


Cix,y) 


-0.23 ±0.11 


-0.26 ±0.04 


0.07 ±0.04 


-0.63 ±0.06 


-0.24 ±0.04 


-0.37 ±0.03 


a 


-o.oitroi 


-0-03t°;f7 


n ono+''oo'' 


_f, -,0+0.06 
'^•^■'-0.10 


-o-o5t°:l^ 


u.04_g4g 


b 


-0.08t°;°l 


r|nfi+0'36 

'J-^0_g2g 


n n"S3+'^oo6 
u.uaci_m)53 


07+''-"2 
'^•'-"-0.01 


n TH+0.20 


70+°-2'' 

'-'■''-'_0.44 


C^int 


o.4ot°:J° 


^••^''-0.05 


,-, qO + 0.05 

U.OO_g Q4 




,-, qo+0.05 

U.OO_g Q4 


■^2+°-°* 



with t;6ar('") the sum of the disk and gas (if any) circular 
velocities. Although such a definition strictly holds only for 
a spherical system and is therefore rigorously wrong for a 
thin disc, we prefer to use it since it can be easily related to 
the observed circular velocity data. 

As for the logarithmic slope, we must choose a reference 
radius to estimate the DM mass fraction. We again consider 
two different possibilities, namely the disc scalelength Rd 
and the halo characteristic radius R-2. We then find: 

ifd) = 0.74 , fd = 0.77 , 68% CL = (0.61,0.89) , 

(/„2> = 0.75 , /-2 = 0.85 , 68% CL = (0.73, 0.91) . 

having defined fd = foAiiRd) and f-2 = foAiiR-i)- Con- 
sistent with the common expectation, it turns out that 
dwarfs galaxies are DM dominated already at the disk scale- 
length radius Rd (and hence still more at R-2 > Rd)- Al- 
though such a result is model dependent, having been ob- 
tained fitting a specific halo model, the agreement of that 
model predictions with the measurements makes us confi- 
dent that the above estimates of {fd, f-2) are fully realistic 
thus confirming the popular picture of dwarfs as ideal labo- 
ratories to study the dark halo properties. 

Since all the dwarfs in our sample are DM dominated 
notwithstanding their different properties, we do not expect 
to find a correlation of {fd, f-2) with {O, Eio, -Dtso)- Table[4] 
shows that both the Spearman correlation coefficient C{x, y) 
and the slope a of the log (fi/fi) vs {Q, log Eio, log O750) are 
strong evidence for the lack of correlation between the DM 
mass content and the environment the dwarf galaxy lives in. 
Actually, according to the 68% CL range for a, one could 
argue in favour of a correlation between fd and the global 
density estimator -D750 , but a look at the data in Fig.[5]shows 
that the points are quite scattered around the best fit line 
and with large errors so that better quality data are needed 
to confirm this preliminary and discrepant result. 



An alternative way to look at the interplay between DM 
halo and baryons is to consider the ratio among a character- 
istic halo size and the disc one. We therefore estimate the 
two quantities 7i—2 = R-2/ Rd (which determines when the 
galaxy starts be fully DM dominated) and TZvir ~ Rvir/Rd 
(which is related to the total mass ratio Mvir/Md)- We find : 

(7^-2) = 5.31 , n-2 = 3.99 , 68% CL = (1.64, 16.75) , 

(7^„ir) = 65.0 , n^ir = 56.2 , 68% CL = (18.8, 166.3) , 

both distributions being characterized by long tails towards 
values larger than the mean ones. Fitting log {TZi/TZi) vs the 
environment estimators, we find the results summarized in 
TableO The slope a of the fitted relations is always well 
consistent with a = within the 68% CL so that we can 
safely conclude that there is no impact of the environment 
on the size ratios. It is also interesting to note that the scat- 
ter in log {TZvir /Tivir) is the same from one case to another 
suggesting that this quantity is roughly constant along the 
sample. Since TZ^ir c>c {A'lvtr/MdY , such a result could be 
anticipated noting that, on such large scale, it is reason- 
able to expect Mmr/Md oc eQ,cDM /^b with Q,cdm {^b) 
the CDM (baryon) cosmological density parameters and e a 
scaling factor quantifying the efficiency of the conversion of 
baryons into stars and HI/Hq gas. 

4.4 Newtonian accelerations 

A quantity that has recently attracted a lot of interest is 
the Newtonian acceleration g{r) = GM{r)/r'^ where M{r) 
is the mass of the DM halo or the baryon component. We 
reconsider here this issue for our dwarfs sample choosing the 
disc scalelength Rd the halo scale radius R-2 as reference 
positions. Setting r — Rd, we find: 



{id > 



-8.55 



-DM 

Id 



-8.48 , 68% CL = (-8.88, -8.22) 
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Table 7. Same as Table[2]but for the Newtonian accelerations at i?_2 using M = 12, 26, 29 galaxies for correlations with {©, Eio, -D750). 



Id 


■yEf-0 


7^2" -Sio 


l^i'-Drso 


7^Y-© 


7^"2'--Sio 


7'"2''-'D750 


C{x,y) 


-0.27 ±0.11 


-0.34 ±0.04 


-0.09 ±0.03 


-0.08 ±0.11 


-0.11 ±0.04 


-0.20 ±0.03 


a 


omt^ol 


-0.15+_IH 


-0.02±°;1? 


nnn+ooi 

U-UU_o.o2 


-0.041°:!? 


1 5+0-24 


b 


0.072tO;«06 


0.24«:^^ 


0.09t«-^ 


o-i8iro2 


-0-04t°i° 


U./tl_g2g 


O'int 


o.26«:l]^ 


0.25«:«^ 


0.251H1 


o.ostH^ 


0.27t«:- 


lO+o"* 

'J-^'J-0.03 



Table 8. Same as Table[2]but for the Newtonian accelerations as a function of (A4B,logi?,;;). 



Id 


I^'-Mb 


l^^'-R^i 


I^/^'-Mb 


7r-«<i 


7?2'^-^B 


7?2"^-i?. 


7^T--^fl 


7^T-^d 


C{x,y) 


-0.09 ±0.04 


-0.19 ±0.04 


0.00 ± 0.04 


-0.53 ±0.03 


-0.08 ±0.04 


-0.27 ±0.03 


-0.39 ±0.03 


-0.13 ±0.04 


a 


"0.02+°:- 


-0 09+°-^^ 


000+"°°'* 


-0 66+'''^* 


-0 01+°-°^ 
'^■'^^-0.02 


-0 24+°-*^ 


-o.o4+°:»t 


-0 02+"'°^ 


b 


-0-55+°-Jg 


-0 04+0-03 

^•'-'^-0.03 


n ns+O-08 


n 14+0-02 


-0 29+°-0^ 


07+°'°^ 

'^-'-"-0.04 


-o.6o+°:« 


no+0-03 

"^■^^-0.03 


O'int 


0.32tro! 


n OO+0.05 


n 90+0.04 


26+°-°'' 


n 90+0.05 
^•"=■^-0.03 


n OO+0.04 


o-o8+°::!l 


15+006 
'^■^'^-0.04 



{7d ) = -9.09 , 7d = -9.10 , 68% CL = (-9.49, -8.70) , 

where we have used the notation 7^ = \og gj{ri) with i = d 
(i = -2) for ri = 7?^ {n = R-2) and j = DM (j = bar) 
for the DM (baryon) Newtonian acceleration (measured in 
cm/s'^). For r — R-2, we get : 



(7-2 > 



-8.62 , 7_2 



-8.62 , 68% CL = (-8.89, -8.29) 



(7-"2'') = -9-37 , 7^T = -9.29 , 68% CL = (-9.74, -9.09) . 

An issue hotly debated concerns whether it is possible to 
define a universal quantity evaluating the DM or baryon ac- 
celeration at a given reference radius. Should this be the 
case, the proposed quantity should take a constant value 
independent on the environment the galaxy lives in. We ex- 
amine this issue by fitting iog[gj{ri) / gj{ri)] as a function 
of (0, logEio,logD75o) for the different ij combinations we 
have considereqj. Tables[S] and [7] summarize the results of 
these fits. For almost all cases, the slope a of the fitted re- 
lation is consistent with a null value within the 68% CL 
indicating the absence of any statistical meaningful corre- 
lation. The only exception is represented by the 7^"*^ which 
seems to correlate with both the tidal index 6 and the global 
density estimator -D750, but surprisingly not with the lo- 
cal density estimator Eiq. The intrinsic scatter is, however, 
quite large and the constraints on a too weak to definitely 
argue in favour of an impact of environment on the baryon 
acceleration at the disc scalelength radius. 

The lack of correlations with the environment and the 
not too large 68% confidence ranges may then be considered 
as evidence that the considered Newtonian accelerations are 
indeed roughly universal quantities. In order to further in- 
vestigate this issue, we perform another series of fits taking 



Note that we define gj{ri) = dex(7^) with dex(a;) = 10^. 



as X variable the absolute B band magnitude Mb and the 
logarithm of the disc scalelength log Rd obtaining the re- 
sults summarized in Table[8l It turns out that the baryon 
acceleration at Rd and the DM one at R-2 significantly cor- 
relates with the disc size so that they can not be considered 
as universal quantity. On the contrary, such a role could be 
played by the DM acceleration at Rd and the Newtonian 
one at R-2 since we have found no correlations with either 
the environment or the galaxy properties. It is worth not- 
ing, however, that the scatter around the median for both 
quantities is not negligible so that one should take this pre- 
liminary conclusion cum grano salts. 



5 CONCLUSIONS 

Being likely dark matter dominated, dwarfs and LSB galax- 
ies are ideal targets to investigate the halo density profile 
in the inner regions. We have here used a large sample of 
such systems with well measured and radially extended ro- 
tation curves to constrain the dark matter halo parameters 
adopting the Einasto model for the density profile. Such a 
model has been proposed as a good description of dark mat- 
ter haloes in high resolution N - body simulations and its via- 
bility has been recently confirmed by its succ ess in fitting ro- 
tatio n curves of the THINGS galaxy sample (jChemin et al. I 
l2011f ). Our result strengthens these findings and allows to 
extend the applicability of the Einasto model to the dwarfs 
galaxies regime which has been not tested by either simula- 
tions or observations. In particular, we find that the Einasto 
slope parameter udm can take values in a wide range which 
overlaps the one predicted by N- body simulations, but with 
a strong preference for shallower profiles. Whether this is 
due to the different mass range probed by dwarfs or a con- 
sequence of the impact of baryon collapse (not taken into 
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account in the simulations) is still an open question which 
deserves further investigation. 

The wide range spanned by udm translates in a simi- 
larly large range for the logarithmic slope a(r) when evalu- 
ated at the disc scalelength radius Rd or at a halo character- 
istic radius G.lR^ir- As we have yet discussed, in order to fit 
the data, the parameters of the Einasto model are adjusted 
in such a way to mimic an almost cored profile in the very 
inner regions and an isothermal or NFW-like in the outer 
regions. How this transition takes place depends on the de- 
tails of the galaxy so that we get a large range for both ad 
and Qo.i)- We can therefore safely conclude that, contrary 
to what was argued by the past N - body simulations, there 
is no observational evidence for any universal halo profile. 

Motivated by this result, one can then wonder what is 
driving the diversity of the halo profiles. In an attempt to 
get some hints on this question, we have here investigated 
the impact of the environment the galaxy lives in on the 
DM properties. To this end, we have resorted to three dif- 
ferent indicators (namely, the tidal index G, the local over- 
density Eio and the global overdensity D750) to quantify 
whether a galaxy is isolated or reside in a rich environment. 
We have then used a Bayesian motivated procedure to fit 
the Einasto parameters and other DM related quantities as 
loglinear functions of (0, Eio, -D750) noting that a null value 
for the slope of the fitted relation is a strong evidence for 
the absence of any correlation. According to this criterium, 
we can argue that the environment does not play any role in 
determining the final properties of the dark matter haloes 
of dwarfs galaxies. Actually, this is somewhat counterintu- 
itive since, e.g., one can easily think of tidal interactions 
playing an important role in modifying the initial density 
profile of the dwarf halo. Should this be the case, we must 
find a difference between the halo parameters according to 
the properties of the environment. Since this is not the case, 
one should look for a different mechanism which is environ- 
ment independent and able to wash out the impact of tidal 
interactions on the dark matter haloes. 

Further work is actually needed in order to strengthen 
the preliminary conclusions we have presented here. First, 
the galaxy sample should be enlarged in order to better con- 
strain the investigated correlations and estimate their intrin- 
sic scatter. Although an updated version of the CNG catalog 
is in preparation (I. Karachentsev, private communication), 
it is actually hard to significantly increase the number of 
galaxies with both an available well measured rotation curve 
and a tidal index determination. As an alternative way, one 
could try to refine the local and global overdensity estima- 
tors in order to reduce the impact of chance projection and 
make them more closely related to the dynamical state of the 
galaxy. Should the combination of an updated sample and 
better environment indicators be available and still confirm 
the absence of any correlation, one should then open a new 
path towards understanding the physical processes driving 
the formation of dark matter haloes and the interactions 
between the dark and luminous galactic components. 



ACKNOWLEDGMENTS 

We warmly thank W.J.G. de Blok, S. Oh, J. Simon and 
R. Swaters for making their data available in electronic for- 



mat. I.D. Karachentsev and V.E. Karachentseva are greatly 
acknowledged for discussion on the use of the tidal index 
and other environment estimators. We are also grateful 
to an anonymous referee for his/her comments that have 
significantly helped to improve the paper. This research 
has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology, under contract with 
the National Areonautics and Space Administration. VFC 
is funded by the Italian Space Agency (ASI). 



APPENDIX A: TESTING THE METHOD 

The analysis presented above relies on the implicit as- 
sumption that our fitting procedure correctly recovers the 
halo model parameters {jidm , c^ir ,^og M^ir) and hence the 
dark matter related quantities of interest we have discussed 
above. Although our method is a standard one, it is still not 
guaranteed a priori that the data probe the rotation curve 
with sufficient accuracy and sampling. Moreover, since the 
data only probe a limited radial range {Rmin , Rmax) , it is 
possible that that Rmin/Rd (Rmax/Rd) is too large (too 
small) to guarantee that we are correctly recovering the in- 
ner slope (the halo mass). It is therefore important to test 
whether our inferred constraints are biased or not. 

To this end, we apply our fitting procedure to a set of 
simulated rotation curves generated according to the steps 
schematically sketched below. 

(i) We choose a reference galaxy from the sample 
used in the text and set the disc {Md,Rd) and gas 
(MhIjRhi) parameters of the simulated galaxy as 
{e^dMd,erdRd,(^mgMHi,ergRHi) with e; randomly gener- 
ated from uniform distributions in the (0.1, 10.1) range for 
the masses and (0.5,2.5) for the scale radii. 

(ii) Concerning the halo model parameters, we first set 
nDM by randomly extracting from a uniform distribution 
over the range (1,11). We then set the virial mass M^ir 
solving the empirically motivated relation : 



Me 



[M,„/(9.5 X 10" Mq)]' 



9 X 109 M0 ; 1 -^ [M„i,/(9.5 X 10" Mq)]-°-^^ ■ 

Finally, we set the concentration Cvir by asking that the halo 
contribution to the total circular velocity at Rd equals udm 
with vdm randomly sampled from the range (0.05, 0.95). It 
is worth stressing that this procedure allows us to simulate 
a wide range of different disc/halo combinations, from disc 
to halo dominated galaxies. Moreover, the halo parameters 
are fully realistic being in agreement with both the result 
of numerical simulations and actual observations. 

(iii) Having set the input baryons and halo parame- 
ters, we sample the theoretical rotation curve over the 
range {£minRd,£maxRd) with 0.1 < £min < 0.6 and 
2.0 ^ Smax ^ 12.0 so that we can mimic different radial 
range coverage. We divide this range in ACjim equally 
spaced bins and take one random point from each bin with 
A/'sim the total number of points in the simulated curve 
chosen to be roughly equal to the one in the reference curve. 
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(iv) We finally assign to each R in the list generated 
above a value of the total circular velocity extracted from 
a Gaussian distribution centred on the theoretical value 
and with a variance set to 0.5%. To this fake observed 
Vc{R), we attach a measurement error which is similar to 
the one on the closest R point in the reference rotation curve. 

As reference galaxy, we use NGC3274, UGC4173 and 
UGC 10310 and generate ~ 50 simulated rotation curves for 
each case thus getting a good statistics to check the efficiency 
of our fitting procedure. For each parameter x, we look at the 
distributions of Ax/a^ with Aa; = Xinp — Xout and {xout, ctx) 
the median and (symmetrized) 68% uncertainty as inferred 
from the fit. Note that, because of the scatter introduced in 
the simulated data, we do not expect to perfectly recover 
the input values (i.e., Aa; = 0). Rather, we can consider the 
fitting procedure as reliable if the Ax/a^ distribution is cen- 
tred on a value smaller than 1. This is indeed the case for all 
the parameters we have considered in the analysis presented 
in the paper. We can briefly summarize the results of this 
test as follows. 

• The halo model parameters (n£)M,c„ir, log M„ir) and 
the stellar M/L ratio are well recovered with the median 
values of Ax/a^ ~ (0.2, 0.3, 0.3, 0.1) and rms values smaller 
than 1 for all of them. Not surprisingly, we find that Ax/a^ 
is anticorrelated with Rmax/Rd for x = (cj,ir,logMjjir), 
i.e., the larger is radial range covered, the better are the 
concentration and virial mass recovered. 

• The logarithmic slope and the dark matter mass 
fraction at the disc scalelength are excellently recovered 
with typical Ax/a^ values close to null for most of the 
cases. Significantly larger values are obtained when the 
same quantities are estimated at the halo radius R-2 as 
a consequence of the uncertainties in recovering both the 
concentration and the virial mass. We, however, stress that 
the median Ax/a^ values are still smaller than unity so 
that we can rely on the estimated values used in the paper. 

• The Newtonian accelerations are correctly recovered 
although with larger A^jo^ values. We have, however, 
checked that there is not any statistically meaningful 
correlation of Ax/gx with any of the parameters entering 
the correlations investigated so that we are confident that 
the estimated quantities are not systematically biased. 

Motivated by these comforting results, we therefore safely 
conclude that our fitting procedure and the available data 
are sufficient to reliably estimate the halo parameters and 
the DM related quantities we have been interested in. 
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